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DERIVATION OF THE DIFFERENTIAL TUNNELING CONDUCTANCE WITHIN THE GUTZWILLER APPROXIMATION 

In this supplemental material, we provide a more detailed discussion about the differential tunneling conductance within the 
Gutzwiller approximation. 

When the STM tip is positioned in the vicinity of site i, the tunneling Hamiltonian is given by: 

Ht = Y,[Vt,{clc^a + H.c.) + Vtf{clf,, + H.c.)] 

(T 

= Y,Vtc{cl[c^. + iVtf/Vtc)h.]+li.C.} . (1) 

(T 

Here cj^. (cta), cf^. (cjcr), fi^ (fia) are the creation (annihilation) operators for the STM tip electrons, and the conduction band 
electron and truly physical /-electrons at site i with spin projection a, while Vtc and Vtf are the amplitudes for tunneling into 
conduction and truly /-electron states. Therefore, the total Hamiltonian consists of three parts: 

-f^total = Htip + i^sample + Ht , (2) 

with Hiip = J2ka ^tkcl^^^ctka represents the Hamiltonian for the STM tip while iJsampie represents the sample system (that is 
the heavy fermion system) as given by Eq. (1) in the main text. The electron tunneling operator is defined as the rate of change 
of the tip electrons fit = electa- 

I = -^["-t, -ff total]- 



^e 



^ Vtc{cl[c,a + {Vtf/Vtc)ha] - H.C.} . (3) 



where e is the electron charge and h is the reduced Planck constant. Within the linear response theory, the expectation value of 
/ is derived as follows ^ : 

+yt,%([i(t),Bt(t')]_)e--^(*"*')-y,,y,^.([it(t),i3(tO]-)e^"^(*-*') 

+yt,l-t/([i3(i),it(iO]->e-*^^(*-*')-ytcV^t/([St(t),i(t')]-)e*'=^(*-*') 
+y,}([B(t),Bt(i')]_)e-^(t-0„y^2([St(i)^S(i')]_)e.ey(t-0} 



h ^ 

= ^{V^^lr^[Ui^{^eV)] + VtcVtslm[U^,^{-eV)] + F*eF*/Im[C/,f.^(-eF)] + Vt)\m[UP,,^ {-eV)]} . (4) 

Here = Y.ka^\<j{t)ctka{t) / ^/l^p 13{t) = Eka /tk<T(Octk^(t)/^^/V^, where A^tip is the normalization factor 

for the tip electron degrees of freedom and c\^{t) — e*^''^™p'°*cj^e~*^'"'p''=*, f}„{t) = e'^'"'p'°*/io-6~'^'"°'"'*' ™d Cta{t) = 
giATiipt^^^g-iK.ipt -j-jjg symbols A'sampie and /<tip denote the Hamiltonian with respect to the respective chemical potentials: 
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Ji'sampie = ^sample ^ MsampieA^sarapie and Ki[p = Fjip - AitipA'tip, wherc iVsampie and iVtip are the total electron numbers in the 
sample and the tip, respectively. The chemical potential difference between the tip and the sample is identified as the voltage 
bias /itip ^ Msampie = . The retarded correlations are defined as; 

U,i■'{t)^~^emA{t),A\0)]_), (5a) 
C/i^(t) = -i0(i)([i(i),Bt(O)]_), (5b) 

Ui''{t) = -ze{t){[B{t),AH0)U) , (5c) 
U,g^{t)^~te{t){[B{t),BH0)U . (5d) 

The Fourier transform of these retarded correlation can be evaluated more easily from their corresponding Matsubara correlation 
functions, which are found to be: 

C/^^(iw„) = X! ^ X] ^iai^Pn)Gt,ka{iPn + i^^m) , (6a) 

C/^^(zc^„0 - X ^ 51 G{^i^Pn)GtMAVn + ioJm) , (6b) 

^ kcr ipn 



U^^{iUJr,i) = -r^ X] ^ X! Gll{ipn)Gt.\^a{iPn + iw^) , (6c) 



f/^^(zw„) = -^Y.\Y. Gf/(^P„)G^,k.(^p„ + zc^„) , (6d) 



where /3 = 1 /ksT with fc^ and T the Boltzmann constant and the electron temperature, respectively, while p„ = (2n7r+ 1)/sbT 
and LUm = 2'rmTkBT are the Matsubara frequency for fermions and bosons, respectively. The electron Green's functions for 
different parts can be written as; 

G(,ko-(«Pn + «w,„) = — , (7 a) 

ipn + lOJm - ?k 

GfA^Pn) = r de^^ , (7b) 

J_oo iPn - e 

GfJ^'\^Pn)^ deP^^^, (7c) 



= d6^^ , (7d) 

J-oo Wn - e 



where we have introduced the local single-particle and hybridization density of states of truly physical electrons: Pi^ie), pfcr(^)' 
and Pii'^^'^\e)- Performing the frequency summation, we obtain the current as: 



AVVtip J 

Within the wide-band approximation for the tip electrons, we arrive at 

(/> = E / d'^^tcPtA^) + mcVtiplUe) + V,)pl{e)][fMe) - h^{e - eV)] . 

where A^o is the density of states for the STM tip. The differential tunneling conductance is then found to be; 

d{J) _ 2e^TTNo ^ f 2 ; 2 f . M ^ g/FD(£-ey)^ 



(8) 



(9) 



. Y: I de[V,lpl{e) + 2V,cV,fptU^) + V,)pl{e)] "^^ ) , (10) 



dV h 
and has the zero-temperature limit 

'^"'"'^ ■ Y.'ytM^V) + 2V,,V,fpfAeV) + V,}pUeV)] . (11) 



d{I) _ 2e^TTNo 



dV 



3 



The same type of formula can be written for a finite temperature when the thermalized local density of states is defined as 
in Eq. (6) of the main text. The remaining task is to evaluate the truly electron local density of states in the sample. It is a 
nontrivial task in strongly correlated electron systems and except for a few rare cases approximations are usually made. In 
the methods like non-crossing approximation (NCA) |0,[3tl, equation-of-motion (EOM) decoupling scheme fllltl, numerical 
renormalization group (NRG) method |@|, or quantum Monte Carlo (QMC) simulation |70, the electron Greens functions are 
directly evaluated. In this fashion, the renormalization from correlation effects is directly encoded into the electron self-energy 
and the above tunneling conductance can be applied directly. However, when we solve the strongly correlated electron systems 
within an auxiliary field theory, a caution must be taken. In the Gutzwiller approximation, we note that the coherent part of 
the Green's functions involving the truly /-electrons are related to the Green's functions involving the quasiparticle operator as 

Gll^uj) = 9iaGll {uj) and (uj) = gf^G{^{uj), where G^^(w) and G{1{uj) are determined by the solutions to the effective 
Hamiltonian given by Eq. (2) in the main text. With these relations, the zero-temperature differential tunneling conductance 
reduces to Eq. (7) in the main text. The same type of formula has also been derived in the large- mean-field approach to the 
Kondo lattice 181]. We conclude that it should be applicable to all theoretical approaches based on an auxiliary field theory (like 
slave-boson mean-field theory |@]). 
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Motivated by recent success of local electron tunneling into heavy fermion materials, we study the local 
electronic structure around a single Kondo hole in an Anderson lattice model and the Fano interference pattern 
relevant to STM experiments. Within the Gutzwiller method, we find that an intragap bound state exists in the 
heavy Fermi liquid regime. The energy position of the intragap bound state is dependent on the on-site potential 
scattering strength in the conduction and /-orbital channels. Within the same method, we derive a new dl /dV 
formulation, which includes explicitly the renormalization effect due to the /-electron correlation. It is found 
that the Fano interference gives asymmetric coherent peaks separated by the hybridization gap. The intragap 
peak structure has a Lorenzian shape, and the corresponding dl/dV intensity depends on the energy location of 
the bound state. 
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PACS numbers; 71.27.+a, 74.55.+V, 75.20.Hr, 75.30.Mb 

Introduction. The intermetallic heavy fermion compounds 
based on either rare earth elements or on actinides |tll43|] ex- 
hibit many unusual properties like heavy Fermi liquid (HFL), 
magnetic ordering, quantum phase transitions and associated 
non-Fermi liquid, as well as unconventional superconduc- 
tivity |4j]. In these materials, there are two types of elec- 
trons: delocalized conduction electrons (c-electrons), which 
derive from outer atomic orbitals, and strongly localized /- 
electrons that singly occupy inner orbitals. It is well ac- 
cepted that the interplay between the c-f hybridization and the 
screened on-site Coulomb repulsion is the key to the above- 
mentioned anomalous properties. However, the precise na- 
ture of the coherent Kondo state responsible for the forma- 
tion of HFLs remains hotly debated i^]^- Recent success of 
scanning tunneling microscopy (STM) experiments on heavy 
fermion materials ||9 - IJ | opens a new avenue toward under- 
standing of these remarkable properties. Interest in the prob- 
lem was also stimulated by the possibility to reveal the nature 
of electronic correlation effects by impurities and defects in 
heavy fermion materials as in unconventional superconduc- 
tors including cuprates 111211 and iron pnictides il3l4l5[l or se- 
lenides lll6ll . 



approach 11 1 8142 111 . However, there is still a significant ques- 
tion about the role of correlation effect in the Fano interfer- 
ence. In one scenario llSll . the passage of an electron from 



the STM tip is accompanied by a simultaneous spin flip of the 
localized moments via cotunneling mechanism. When the lo- 
cal spin operator is represented in terms of pseudofermions, 
this cotunneUng term is equivalent to an effective tunneling 
in the pseudofermion channel renormalized by a local boson 
condensation order parameter inherent to the Kondo lattice 
itself ilSll . In another scenario within the same Kondo lat- 
tice model |fl9!l, the renormalization factor is absent in the 
dl/dV formula. For both the single Kondo impurity problem 
and the Kondo lattice model, since the renormalization factor 
can be absorbed into the bare tunneling amplitude, this differ- 
ence will not cause qualitative difference in analysis of STM 
data. However, for a Kondo hole or a more general disordered 
Kondo lattice problem, a logically consistent dl/dV formula- 
tion is important. In connection with the STM measurements, 
the latter type of problems just began to attract more focused 
interest because the inhomogeneity may provide important in- 
sight into the complex electronic structure of heavy fermion 
materials. 



Theoretical challenge in the interpretation of differential 
tunneling conductance as measured by STM in heavy-fermion 
lattice systems and around a single magnetic impurities or, 
more generally, disordered Kondo lattice systems, requires 
a proper treatment of electron correlation effects and quan- 
tum interference between the electrons tunneling from the 
STM tip into the conduction band and into the magnetic /- 
electron states. In the single Kondo impurity case, the line 
shape of dl/dV has been described reasonably well with a 
phenomenological form, as first discussed by Fano ||17|] . A 
microscopic understanding of this Fano line shape was re- 
cently provided | ]_8, 3 in the single Kondo impurity case. 
In addition, the essence of a similar line shape observed on 
a clean crystal surface which has the lattice transla- 

tional symmetry, can be captured within a similar microscopic 



The Kondo hole problem was previously studied within an 
Anderson lattice model, where the /-electron self-energy for 
the pristine system was obtained from the second-order per- 
turbation in Hubbard repulsion ll22ll . In the Kondo insula- 
tor (KI) regime, in which a hybridization gap is open at the 
Fermi energy, it was shown from the /-electron local den- 



sity of states II22II that the Kondo hole introduced an intragap 
bound state. However, a very recent study of a single Kondo 
hole problem within the Kondo lattice model indicated that 
the existence of the intragap bound state depends on the na- 
ture of defects in the HFL regime 12311 . This discrepancy sug- 
gests further theoretical studies are needed. In addition, exist- 
ing studies I221 12311 are limited to local density of states and 
an analysis of Fano interference in the Kondo hole problem 
is still lacking. In this Letter, we study the local electronic 
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structure around a single Kondo hole in an Anderson lattice 
model as well as the Fano interference pattern relevant to STM 
experiments. Within the Gutzwiller method, we are able to 
demonstrate that the intragap bound state exists regardless of 
whether the system is in the HFL or KI regime. The energy 
position of the intragap bound state is dependent on the on-site 
potential scattering strength in the conduction and /-orbital 
channels. The same Gutzwiller method enables us to include 
explicitly the renormalization effect due to the /-electron cor- 
relation into the dl/dV formulation, which is conceptually 
consistent with the cotunneling mechanism 1 1 80 . It is found 
that the Fano interference does give rise to an asymmetric co- 
herent peaks separated by the hybridization gap, while it tunes 
the intragap peak structure dependent on the energy location 
of the bound state, as shown in dl/dV. 

Kondo hole model and Gutzwiller formalism. Our Hamil- 
tonian for the single Kondo hole system is written as 



E 



H.c. 



(1) 



Here the operators cj^ (cia) create (annihilate) a conduction 
electron at site with spin projection a while the operators 
firr ifia) Create (annihilate) a /-level electron at site with 
spin projection a. The number operators for c and / orbitals 
with spin a are given by n,^ = 4^Cia and = f^^fta, re- 
spectively. The quantity f^^ is the hopping integral of the con- 
duction electrons, and e/ is the local /-orbital energy level 
on the magnetic atoms. The hybridization between the con- 
duction and /-orbital on the magnetic atoms is represented 
by Vcf and the /-electrons on the magnetic atoms experience 
the Coulomb repulsion of strength U. The first two lines on 
the right-hand side of Eq. (|T|) constitute the standard Ander- 
son lattice model; while the last two lines represent the effect 
of a doped nonmagnetic atom, which without loss of general- 
ity is located at the origin vj = (0,0). The /-orbital energy 
level on the singly doped nonmagnetic atom is given by In 
the Kondo hole problem, ej^ will be adjusted to ensure there 
is no /-electron occupation on the missing /-character center 
and as such the effect of on-site Coulomb repulsion can be 
negligible in the study of low energy electronic structure. In 
addition, the Kondo hole will also give rise to a potential scat- 
tering potential and a possible change of local hybridization 
Vj^ . These three parameters make the description of a Kondo 
hole more realistic, demonstrating the flexibility of the An- 
derson lattice Hamiltonian. For simplicity, the effect of local 
hybridization change is neglected by assuming Vjj — Vcf in 
the present work. 

Due to the presence of onsite Hubbard interaction U be- 
tween the /-electrons on the magnetic atoms in Eq. the 



above problem is strongly correlated. This strong correlation 
effect can be accounted for by reducing the statistical weight 
of double occupation in the Gutzwiller projected wavefunc- 
tion approach 12411 . and the projection can be carried out semi- 
analytically within the Gutzwiller approximation l25l - l27ll . In 
the present problem, the lattice translation symmetry is bro- 
ken due to the Kondo hole, we use a spatially unrestricted 
Gutzwiller approximation (SUGA) l28l - l32ll to translate the 
original Hamiltonian Eq. ([T]) into the following renormalized 
mean-field Hamiltonian: 

ij,(T iT^-^-cr 

+ E + - ^')flf^- + E + E ^ccLc/. 



+ E(4 - + Y}^lAJi^ + H.c] , (2) 

where A^cr and di are the Lagrange multiplier and the double 
occupation at site i. We have used fj^ {fia) to denote the 
quasiparticle field operators to differentiate from the truly /- 
electron operators in Eq. ([T]). The local c-f hybridization has 
been renormalized by a factor of gi^, which is given by 



{nl - di){l ~ nl + di) 



f 



1/2 



1/2 



(3) 

with being the expectation value of the spin-tr density 



operator = fj^fia and hj ~ X^cr^L- Minimization 
of the expectation value of Hgff leads to the following self- 
consistency conditions for A^cr and df. 



-U 



^c/E%^((4/..)+c.c) 



dd,. 



(4a) 
(4b) 



for i ^ 1. Equation (|2]i can be cast into the Anderson- 
Bogoliubov-de Gennes (Anderson-BdG) equations 133 1 




En 



(5) 



subject to the constraints given by Eq. (HJi. Here /i^^ = 

and h{^ = [(e/ + Xi){l - dii) + ejSu - fi]Sij. Since most 
of heavy-fermion systems have a layered structure, in which 
the dominant effects occur in the planes containing /-electron 
atoms, we solve this set of equations self-consistently via ex- 
act diagonalization on a two-dimensional square lattice. After 
the self-consistency is achieved, one can then calculate the 
projected local density of states (LDOS) as defined by: 



dE 



(6) 
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FIG. 1. (Color online) Solution to the Kondo hole problem in the 
heavy Fermi liquid regime for ^ — —0.4. Contour plots of the 
Gutzwiller factor gi (a), the Lagrange multiplier Xi, the double oc- 
cupation di (c), the c-electron density n1{d), /-electron density n{ 

(e), and the hybridization density n^-'^ (f). The parameter values 
e'f — 100 and = 0. The values of other parameters are given 
in the main text. 

where the Fermi-Dirac distribution function fpoiE) = 
[cxp{E / ksT) + Here we have used the fact that the 

eigenfunctions are real and the system has a two-fold spin de- 
generacy in the non-magnetic state, for which quantities like 

n{^, gia, and Xi„ become spin-independent. Throughout this 
work, the quasiparticle energy is measured with respect to the 
Fermi energy and the energy unit = 1 is chosen. 

Local electronic structure around the Kondo hole. In our 
numerical calculations, we take the following values of pa- 
rameters for the pristine system: The on-site Coulomb inter- 
action on /-electrons U = 2, the bare hybridization Vcf = 1, 
and the local /-level is taken to be e/ = —1. The temperature 
is fixed at T = 0.01. We solve the Anderson-BdG equations 
on a 32 X 32 square lattice (by assuming a periodic bound- 
ary condition) to determine the self-consistency parameters gi 
and \i while calculate the LDOS by using the supercell tech- 
nique |l34[| (8x8 supercells are used). The chemical potential 
is varied so that the pristine system can be tuned into the HFL 
or KI regime. To describe the Kondo hole, we have fixed the 
value of to be positive. This choice is consistent with the 
experimental realization of Kondo holes by substituting a Ce 
ion in a stoichiometric Ce compound by a La atom or a U ion 
in a U-heavy-fermion system by a Th atom, where the /-level 
is unoccupied by electrons on these impurity atoms. 

Figure [T] shows the spatial variation of the self-consistently 



determined Gutzwiller factor gi (panel (a)), Lagrange multi- 
plier \i (panel (b)), the double occupation di, as well as partial 

— — f — cf 

charge density n"^ (panel (d)), ni (panel (e)), and n,- (panel 
(f)) around the Kondo hole in the HFL regime with fi — —0.4. 
Here = 100 and = 0. The change of all these quantities 
happens around the Kondo hole. The short length scale of the 
change is consistent with the coherence length as estimated 
by ^ = h{vp) /ngVcf, where (vp) and g are the averaged 
Fermi velocity and the Gutzwiller renormalization factor for 
the pristine system. With the given set of parameter values, we 
estimate h{vp) = 4.89, and g ~ 0.95 from the self-consistent 
iterations, which leads to ^ being about 1.64 lattice constants. 
The anisotropy of the spatial change follows the underlying 
lattice, which is similar to that expected for a superconduct- 
ing order parameter around a unitary non-magnetic impurity 
in a short-coherence s-wave superconductor With the cho- 
sen parameter, it is found that the Kondo hole is negatively 
charged in the conduction band, by which we mean the con- 
duction electron density on the Kondo hole is smaller than the 
value for the underlying pristine system. To get a positively 
charged Kondo hole in the conduction band, an attractive po- 
tential scattering must be introduced. 

In Fig. |2] we show the local partial and hybrid DOS in the 
HFL regime but with various values of and [The profile 
of self-consistency quantities around the Kondo hole remains 
similar to those shown in Fig. [T] and not shown here]. In the 
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FIG. 2. (Color online) Local density of states, in the HFL regime, on 
the Kondo hole site (left column) and on its nearest neighboring site 
(right column). Panels with label (aX), (bX), and (cX) correspond 
to the conduction band and /-quasiparticle, and hybrid density of 
states, respectively. The red-dashed line is for the pristine case; the 
black-solid and blue-dotted lines for €f — 1.0 and 100 with el — 
fixed, the green-dashdotted line for ey = 100 while = —1.0. The 
coherent peaks for the pristine system are located at i5 = 0.016 and 
0.504. The bound state peak due to the presence of Kondo hole is 
located at S = -0.082 and 0.224 for = 1 and 100 with = 
fixed, and at E = 0.064 for ei = 100 but ei = -1.0. 
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FIG. 3. (Color online) dl /dV, in the heavy Fermi liquid regime, as a 
function of energy on the Kondo hole (a) and its nearest neighboring 
site (b) corresponding to Fig. (2] as well as the dl/dV spatial map 
with selected energies E = 0.016 (c), 0.224 (d), and 0.504 (e). These 
energies correspond to the coherent peak position (E = 0.016 and 
0.504) and the intra-gap peak position (E = 0.224) for ei 



100 



and e„ 



0. 



absence of Kondo hole, the LDOS is spatially independent, 
but both the local partial and hybrid DOS is asymmetric with 
respect to the chemical potential (E = 0). The DOS intensity 
is peaked at the coherent gap edge (see the red-dashed line in 
the figure) and is finite at the Fermi energy, indicating the elec- 
tron states are in the HFL regime. In the presence of Kondo 
hole, when ej^ is increased, an intragap bound state is formed 
with the energy position moved toward the center of the gap. 
With increased the DOS intensity on the /-channel van- 
ishes directly on the Kondo hole site while becomes stronger 
on the site nearest neighboring to the Kondo hole. However, 
the corresponding DOS in the conduction band has an oppo- 
site trend. Furthermore, an attractive potential scattering on 
the conduction channel will further tune the bound state peak 
away from the gap center It is also interesting to note that 
the hybrid DOS exhibits negative-sign peak for the intragap 
state for weak to intermediate values of ej^ and e^. Our result 
agrees with the STM observation of an electronic bound state 
at thorium atoms, when they are substituted for U atoms in 
URusSia Q. 

Fano interference pattern around the Kondo hole. To 
model the STM electron tunneling, for which a schematic 
picture can be found in Fig. 1 of Ref. we introduce the 
tunneling Hamiltonian between the STM tip and the sample 
at a specified measure site "z": Ht ~ X^cr ^tc{cta[cicr + 
{Vtf /Vtc),fia] + H.c.}, where Vtc and Vtf are the amplitudes 
for tunneling into conduction and /-electron states. Within 
the Gutzwiller approach to account for the /-electron corre- 
lation, tunneling amplitude Vt / will be renormalized by the 
local Gutzwiller factor gi, that is, Vtf ^ diVtf, which leads 
to, in the weak-tunneling limit, the differential tunneling con- 
ductance 13611: 



where r = Vtf /Vtc- Here we have assumed the tip density of 



states to be independent of energy and approximated its value 
at the Fermi energy Nq, which is on the order of magnitude of 
the inverse band width of the tip. Equation (|7]i indicates that 
the renormalization due to the /-electron correlation effect in 
the heavy fermion materials must be taken into account ex- 
plicitly in the tunneling conductance formula. This formula is 
logically consistent with the cotunneling mechanism llST . and 
must be used in theoretical approaches based on an auxiliary 
field theory. 

As is shown in Fig. |2] the hybrid DOS is not positively 
definite, the line shape of dl / dV is determined not only by 
whether the electronic structure for the pristine system has 
the particle-hole symmetry but also by the relative strength of 
tunneling amplitude of conduction band and /-electron chan- 
nel. Figure [3] shows the energy dependence of dl/dV on the 
Kondo hole (panel (a)) and its nearest-neighboring site (panel 
(b)), as well as its spatial dependence at selected energies 
(panels (c)-(e)), in the HFL regime with r = 0.2. We used 
the conductance formula (Eq. for T = 0.01 as a very low 
temperature solely for technical reasons. As can be seen from 
Fig. [3ta)-(b), the Fano interference makes the dl/dV char- 
acteristic strongly asymmetric, with the the continuum part 
of dl/dV intensity at negative energies much smaller than at 
positive energies. The Fano interference also leads to different 
dl/dV intensity maps at the two gap edges (see Fig.[3|c) and 
(e)) with a ripple-like structure easily visible at the positive 
energy side. However, the peak structure due to the intragap 
bound state formed around the Kondo hole is robust against 
the Fano interference. On the Kondo hole site, the overall 
dl/dV characteristic is similar to that of LDOS on the con- 
duction band; while its nearest neighboring site, the charac- 
teristic is sensitive to the detailed parameter values of and 
e^. In addition, the spatial imaging of dl/dV characteristic 
(Fig- Bd)) suggests the Kondo hole induced intragap states 
are localized states. 

Conclusion. We have studied the local electronic structure 
around a single Kondo hole in an Anderson lattice model and 
the Fano interference pattern relevant to STM experiments. 
Within the Gutzwiller method, we have obtained the exis- 
tence of the intragap bound state induced around the Kondo 
hole in the HFL regime. The energy position of the intra- 
gap bound state is dependent on the on-site potential scatter- 
ing strength in the conduction and /-orbital channels. Within 
the same method, we have also derived a new dl/dV for- 
mulation, which includes explicitly the renormalization effect 
due to the /-electron correlation. It has been found that the 
Fano interference gives rise to highly asymmetric coherent 
peak structure separated by the hybridization gap. The intra- 
gap peak structure has a Lorenzian shape. The corresponding 
dl/dV intensity depends on the energy location of the bound 
state and the location of the STM measuring point. We have 
also addressed the Kondo hole problem in the Kondo insula- 
tor regime with /i = (not shown). The major finding for the 
HFL regime still holds except for the fact that the electronic 
states are fully gapped at the Fermi energy, and the Kondo 
hole becomes positively charged for = 0. 
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